The Cloning and Characterization of a Three-Finger Toxin Homolog (NXH8) from the Coralsnake Micrurus corallinus That Interacts with Skeletal Muscle Nicotinic Acetylcholine Receptors.

Coralsnakes (Micrurus spp.) are the only elapids found throughout the Americas. They are recognized for their highly neurotoxic venom, which is comprised of a wide variety of toxins, including the stable, low-mass toxins known as three-finger toxins (3FTx). Due to difficulties in venom extraction and availability, research on coralsnake venoms is still very limited when compared to that of other Elapidae snakes like cobras, kraits, and mambas. In this study, two previously described 3FTx from the venom of M. corallinus, NXH1 (3SOC1_MICCO), and NXH8 (3NO48_MICCO) were characterized. Using in silico, in vitro, and ex vivo experiments, the biological activities of these toxins were predicted and evaluated. The results showed that only NXH8 was capable of binding to skeletal muscle cells and modulating the activity of nAChRs in nerve-diaphragm preparations. These effects were antagonized by anti-rNXH8 or antielapidic sera. Sequence analysis revealed that the NXH1 toxin possesses eight cysteine residues and four disulfide bonds, while the NXH8 toxin has a primary structure similar to that of non-conventional 3FTx, with an additional disulfide bond on the first loop. These findings add more information related to the structural diversity present within the 3FTx class, while expanding our understanding of the mechanisms of the toxicity of this coralsnake venom and opening new perspectives for developing more effective therapeutic interventions.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes.

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes.

BACKGROUND: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). RESULTS: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. CONCLUSIONS: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region.

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

miRNAs derived from cobra venom exosomes contribute to the cobra envenomation.

Currently, there is an increasing amount of evidence indicating that exosomes and the miRNAs they contain are crucial players in various biological processes. However, the role of exosomes and miRNAs in snake venom during the envenomation process remains largely unknown. In this study, fresh venom from Naja atra of different ages (2-month-old, 1-year-old, and 5-year-old) was collected, and exosomes were isolated through ultracentrifugation. The study found that exosomes with inactivated proteins and enzymes can still cause symptoms similar to cobra envenomation, indicating that substances other than proteins and enzymes in exosomes may also play an essential role in cobra envenomation. Furthermore, the expression profiles of isolated exosome miRNAs were analyzed. The study showed that a large number of miRNAs were co-expressed and abundant in cobra venom exosomes (CV-exosomes) of different ages, including miR-2904, which had high expression abundance and specific sequences. The specific miR-2094 derived from CV-exosomes (CV-exo-miR-2904) was overexpressed both in vitro and in vivo. As a result, CV-exo-miR-2904 induced symptoms similar to cobra envenomation in mice and caused liver damage, demonstrating that it plays a crucial role in cobra envenomation. These results reveal that CV-exosomes and the miRNAs they contain play a significant regulatory role in cobra envenomation. Our findings provide new insights for the treatment of cobra bites and the development of snake venom-based medicines.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily.

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Comparative Venom Multiomics Reveal the Molecular Mechanisms Driving Adaptation to Diverse Predator-Prey Ecosystems in Closely Related Sea Snakes.

Predator-prey arms races are ideal models for studying the natural selection and adaptive evolution that drive the formation of biological diversity. For venomous snakes, venom is a key bridge linking snakes with their prey, but whether and how venom evolves under the selection of diet remains unclear. Here, we focused on two closely related sea snakes, Hydrophis cyanocinctus and Hydrophis curtus, which show significant differences in prey preferences. Data-independent acquisition (DIA)-based proteomic analysis revealed different degrees of homogeneity in the venom composition of the two snakes, which was consistent with the differential phylogenetic diversity of their prey. By investigating the sequences and structures of three-finger toxins (3FTx), a predominant toxin family in elapid venom, we identified significant differences between the two sea snakes in the binding activity of 3FTx to receptors from different prey populations, which could explain the trophic specialization of H. cyanocinctus. Furthermore, we performed integrated multiomic profiling of the transcriptomes, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and proteomes of the venom glands; constructed venom-related mRNA-miRNA-lncRNA networks; and identified a series of noncoding RNAs involved in the regulation of toxin gene expression in the two species. These findings are highly informative for elucidating the molecular basis and regulatory mechanisms that account for discrepant venom evolution in response to divergent diets in closely related snakes, providing valuable evidence for the study of coselection and coevolution in predator-prey ecosystems.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.).

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Polymerase chain reaction-based snake origin tracing in commercial venom crystals by targeting the mitochondrial D-loop.

Snake venom is a valuable raw material for numerous therapeutic formulations because of its life-saving pharmacological potential. However, due to their high price, fake "snake venoms" have captured a significant portion of the global market, and there is currently no reliable reported DNA-based method available for quickly distinguishing between fakes and originals. Therefore, in this study, a set of newly designed snake-specific universal primers targeting mitochondrial D-loop fragments were employed to detect snake origins in commercial venom crystals by only simplex polymerase chain reaction analysis. Under the optimal thermal cycling conditions, only the 145-149 bp snake-specific mitochondrial D-loop fragments from pure and mixed backgrounds were amplified by the newly designed primers. Specificity was achieved by confirming no DNA amplification occurred in the DNA admixture of ten different chordates, and universality by individual DNA amplification of nine different snakes. The primers that efficiently amplified the minimum mitochondrial DNA contained in a total of 10-2 ng in a 10.0 µl reaction were also successfully able to detect the snake origin in commercial cobra venom crystals. These findings suggest that the newly designed primers can be used to differentiate the original and fake commercial snake venom crystals in order to achieve the highest standards of snake venom-based medications through amplifying the snake-specific mitochondrial D-loop fragments.

Virus-like particles displaying conserved toxin epitopes stimulate polyspecific, murine antibody responses capable of snake venom recognition.

Antivenom is currently the first-choice treatment for snakebite envenoming. However, only a low proportion of antivenom immunoglobulins are specific to venom toxins, resulting in poor dose efficacy and potency. We sought to investigate whether linear venom epitopes displayed on virus like particles can stimulate an antibody response capable of recognising venom toxins from diverse medically important species. Bioinformatically-designed epitopes, corresponding to predicted conserved regions of group I phospholipase A2 and three finger toxins, were engineered for display on the surface of hepatitis B core antigen virus like particles and used to immunise female CD1 mice over a 14 weeks. Antibody responses to all venom epitope virus like particles were detectable by ELISA by the end of the immunisation period, although total antibody and epitope specific antibody titres were variable against the different epitope immunogens. Immunoblots using pooled sera demonstrated recognition of various venom components in a diverse panel of six elapid venoms, representing three continents and four genera. Insufficient antibody yields precluded a thorough assessment of the neutralising ability of the generated antibodies, however we were able to test polyclonal anti-PLA2 IgG from three animals against the PLA2 activity of Naja nigricollis venom, all of which showed no neutralising ability. This study demonstrates proof-of-principle that virus like particles engineered to display conserved toxin linear epitopes can elicit specific antibody responses in mice which are able to recognise a geographically broad range of elapid venoms.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins.

BACKGROUND: The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes. RESULTS: Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae - Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae - Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms. CONCLUSIONS: We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator-prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Cytotoxin antibody-based colourimetric sensor for field-level differential detection of elapid among big four snake venom.

Development of a rapid, on-site detection tool for snakebite is highly sought after, owing to its clinically and forensically relevant medicolegal significance. Polyvalent antivenom therapy in the management of such envenomation cases is finite due to its poor venom neutralization capabilities as well as diagnostic ramifications manifested as untoward immunological reactions. For precise molecular diagnosis of elapid venoms of the big four snakes, we have developed a lateral flow kit using a monoclonal antibody (AB1; IgG1 - κ chain; Kd: 31 nM) generated against recombinant cytotoxin-7 (rCTX-7; 7.7 kDa) protein of the elapid venom. The monoclonal antibody specifically detected the venoms of Naja naja (p < 0.0001) and Bungarus caeruleus (p<0.0001), without showing any immunoreactivity against the viperidae snakes in big four venomous snakes. The kit developed attained the limit of quantitation of 170 pg/µL and 2.1 ng/µL in spiked buffer samples and 28.7 ng/µL and 110 ng/µL in spiked serum samples for detection of N. naja and B. caeruleus venoms, respectively. This kit holds enormous potential in identification of elapid venom of the big four snakes for effective prognosis of an envenomation; as per the existing medical guidelines.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon.

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 µg/g) and long-chain neurotoxin (i.p. 0.14 µg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization.

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

Two Reference-Quality Sea Snake Genomes Reveal Their Divergent Evolution of Adaptive Traits and Venom Systems.

True sea snakes (Hydrophiini) are among the last and most successful clades of vertebrates that show secondary marine adaptation, exhibiting diverse phenotypic traits and lethal venom systems. To better understand their evolution, we generated the first chromosome-level genomes of two representative Hydrophiini snakes, Hydrophis cyanocinctus and H. curtus. Through comparative genomics we identified a great expansion of the underwater olfaction-related V2R gene family, consisting of more than 1,000 copies in both snakes. A series of chromosome rearrangements and genomic structural variations were recognized, including large inversions longer than 30 megabase (Mb) on sex chromosomes which potentially affect key functional genes associated with differentiated phenotypes between the two species. By integrating multiomics we found a significant loss of the major weapon for elapid predation, three-finger toxin genes, which displayed a dosage effect in H. curtus. These genetic changes may imply mechanisms that drove the divergent evolution of adaptive traits including prey preferences between the two closely related snakes. Our reference-quality sea snake genomes also enrich the repositories for addressing important issues on the evolution of marine tetrapods, and provide a resource for discovering marine-derived biological products.

The omega-loop of cobra cytotoxins tolerates multiple amino acid substitutions.

Cobra cytotoxins (CTs), the three-fingered proteins, feature high amino acid sequence homology in the beta-strands and variations in the loop regions. We selected a pair of cytotoxins from Naja kaouthia crude venom to clarify the sequence-structure relationships. Using chromatography and mass spectroscopy, we separated and identified the mixture of cytotoxins 2 and 3, differentiated by the only Val 41/Ala 41 substitution. Here, using natural abundance 13C, 15N NMR-spectroscopy we performed chemical shift assignments of the signals of the both toxins in aqueous solution in the major and minor forms. Combining NOE and chemical shift data, the toxins' spatial structure was determined. Finally, we proved that the tip of the "finger"-2, or the loop-2 of cytotoxins adopts the shape of an omega-loop with a tightly-bound water molecule in its cavity. Comparison with other NMR and X-ray structures of cytotoxins possessing different amino acid sequences reveals spatial similarity in this family of proteins, including the loop-2 region, previously considered to be flexible.

De Novo Venom-Gland Transcriptomics of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia-Next-Generation Sequencing, Functional Annotation and Toxinological Correlation.

Envenomation resulted from sea snake bite is a highly lethal health hazard in Southeast Asia. Although commonly caused by sea snakes of Hydrophiinae, each species is evolutionarily distinct and thus, unveiling the toxin gene diversity within individual species is important. Applying next-generation sequencing, this study investigated the venom-gland transcriptome of Hydrophis curtus (spine-bellied sea snake) from Penang, West Malaysia. The transcriptome was de novo assembled, followed by gene annotation and sequence analyses. Transcripts with toxin annotation were only 96 in number but highly expressed, constituting 48.18% of total FPKM in the overall transcriptome. Of the 21 toxin families, three-finger toxins (3FTX) were the most abundantly expressed and functionally diverse, followed by phospholipases A2. Lh_FTX001 (short neurotoxin) and Lh_FTX013 (long neurotoxin) were the most dominant 3FTXs expressed, consistent with the pathophysiology of envenomation. Lh_FTX001 and Lh_FTX013 were variable in amino acid compositions and predicted epitopes, while Lh_FTX001 showed high sequence similarity with the short neurotoxin from Hydrophis schistosus, supporting cross-neutralization effect of Sea Snake Antivenom. Other toxins of low gene expression, for example, snake venom metalloproteinases and L-amino acid oxidases not commonly studied in sea snake venom were also identified, enriching the knowledgebase of sea snake toxins for future study.

Electric Blue: Molecular Evolution of Three-Finger Toxins in the Long-Glanded Coral Snake Species Calliophis bivirgatus.

The genus Calliophis is the most basal branch of the family Elapidae and several species in it have developed highly elongated venom glands. Recent research has shown that C. bivirgatus has evolved a seemingly unique toxin (calliotoxin) that produces spastic paralysis in their prey by acting on the voltage-gated sodium (NaV) channels. We assembled a transcriptome from C. bivirgatus to investigate the molecular characteristics of these toxins and the venom as a whole. We find strong confirmation that this genus produces the classic elapid eight-cysteine three-finger toxins, that δδ-elapitoxins (toxins that resemble calliotoxin) are responsible for a substantial portion of the venom composition, and that these toxins form a distinct clade within a larger, more diverse clade of C. bivirgatus three-finger toxins. This broader clade of C. bivirgatus toxins also contains the previously named maticotoxins and is somewhat closely related to cytotoxins from other elapids. However, the toxins from this clade that have been characterized are not themselves cytotoxic. No other toxins show clear relationships to toxins of known function from other species.

Mambalgin-1 pain-relieving peptide locks the hinge between α4 and α5 helices to inhibit rat acid-sensing ion channel 1a.

Acid-sensing ion channels (ASICs) are proton-gated cationic channels involved in pain and other processes, underscoring the potential therapeutic value of specific inhibitors such as the three-finger toxin mambalgin-1 (Mamb-1) from snake venom. A low-resolution structure of the human-ASIC1a/Mamb-1 complex obtained by cryo-electron microscopy has been recently reported, implementing the structure of the chicken-ASIC1/Mamb-1 complex previously published. Here we combine structure-activity relationship of both the rat ASIC1a channel and the Mamb-1 toxin with a molecular dynamics simulation to obtain a detailed picture at the level of side-chain interactions of the binding of Mamb-1 on rat ASIC1a channels and of its inhibition mechanism. Fingers I and II of Mamb-1 but not the core of the toxin are required for interaction with the thumb domain of ASIC1a, and Lys-8 of finger I potentially interacts with Tyr-358 in the thumb domain. Mamb-1 does not interfere directly with the pH sensor as previously suggested, but locks by several contacts a key hinge between α4 and α5 helices in the thumb domain of ASIC1a to prevent channel opening. Our results provide an improved model of inhibition of mammalian ASIC1a channels by Mamb-1 and clues for further development of optimized ASIC blockers.

Convergent evolution of pain-inducing defensive venom components in spitting cobras.

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.